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Abstract. We present next-to- leading order QCD corrections to Higgs production in 
association with a photon via weak boson fusion at a hadron collider. Utilizing the fully 
flexible parton level Monte-Carlo program VBFNLO, we find small overall corrections, while the 
shape of some distributions is sensitive to radiative contributions in certain regions of phase- 
space. Residual scale uncertainties at next-to-leading order are at the few-percent level. Being 
perturbatively well under control and exhibiting kinematic features that allow to distinguish it 
from potential backgrounds, this process can serve as a valuable source of information on the 
Hbb Yukawa coupling. 



1. Introduction 

Even more than forty years after the incorporation of the Higgs mechanism into the theory of 
electroweak interactions by Glashow and Salam, its key ingredient, the Higgs boson, still awaits 
discovery. In this context, physicists around the world hope for the Large Hadron Collider 
(LHC) to shed light on the nature of symmetry breaking. 

A particularly promising class of reactions in the context of the search for the Higgs boson 
is provided by weak vector boson fusion processes (WBF). Higgs production via WBF mainly 
proceeds via the scattering of two quarks by the exchange of massive gauge bosons in the t- 
channel which in turn produce a Higgs boson. While the decay products of the Higgs boson 
are located in the central-rapidity range of the detector, the scattered quarks give rise to two 
jets in the far forward and backward regimes, thereby providing a very distinct event topology. 
Typically, these tagging jets are widely separated in rapidity and, due to the electroweak nature 
of the gauge boson exchange, additional QCD radiation is emitted preferentially close to one of 
the original partons in the hard process. This is different from the characteristics of the major 
QCD backgrounds, whose colour structur leads to widespread hadronic activity [1]. 

If the Higgs boson exists, WBF will not only serve as a potential discovery channel, but also 
allow to measure its mass, CP properties and couplings to the top quark, tau lepton and weak 
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gauge bosons [21 13 HI [5]. However, the determination of the Hbb Yukawa couphng remains 
difficult [6j. The H ^ bb decay channel is overwhelmed by large QCD backgrounds which 
are dominated by gluonic components (Fig. [T]). New strategies were proposed to cope with 
these difficulties, such as exploring the jet substructure in so-called "fat jets" [7j. Alternatively, 
requiring the presence of a boson in association with the Higgs boson in WBF processes can 
help to constrain the bottom quark Yukawa coupling [SJ E]. The major drawback of the extra 
heavy gauge boson emission is a significant decrease of the cross section. This loss in statistics 
can be avoided to some extent by requiring the radiation of a photon instead of a boson [lOj. 

Most importantly, the extra photon requirement gives rise to a significant suppression of the 
gluon-dominated backgrounds, due to the absence of a gluon-to-photon coupling in the context 
of the Standard Model. Furthermore, in quark-scattering contributions that are connected by 
a neutral t-channel gauge boson, destructive interference occurs between diagrams emitting a 
photon from an incoming or an outgoing fermion line, respectively. This effect arises in gluon- 
mediated QCD backgrounds as well as in neutral current contributions to WBF processes, such 
as the Z-exchange graphs of Fig. [2] (b). Therefore, the requirement of an additional hard central 
photon does not only suppress QCD backgrounds, but also the ZZ fusion component of the 
signal channel, while the charged-current WBF contributions are relatively enhanced. 

A detailed analysis of the H'jjj signal and its major backgrounds has been performed based 
on leading-order (LO) simulations in Ref. [10]. In that work it has been shown that the central 
photon requirement reduces the cross section of the Higgs signal in WBF approximately by 
a factor of 1/100, as expected from a naive estimate based on the size of the electromagnetic 
coupling constant. The relevant backgrounds instead drop by a factor of about 1/3000, resulting 
in a statistical significance of 1 ^ S/^/B ^ 3 for a luminosity of 100/fb at a center-of-mass energy 
of ^/s = 14 TeV and a Higgs mass of rriH = 120 GeV. In a more recent parton shower study [1] 
a significance of the same order of magnitude has been obtained. 

However, in order to obtain precise predictions accompanied by reliable estimates of their 
theoretical uncertainties, a next-to-leading order (NLO) QCD calculation of the signal process 
is essential. In [11] . we have provided such a calculation within the parton level Monte-Carlo 
program VBFNLO |12^ 113] . which is designed for the simulation of processes with electroweak 
bosons at hadron colliders. 




2. Details of the calculation 

2.1. Tree-level calculation and approximations 

WBF-induced Higgs boson production in association with a photon at hadron colhders mainly 
proceeds via quark-quark scattering processes, qq — )• qqH'j, mediated by the exchange of a 
massive weak boson in the t-channel. The Higgs boson is radiated off this weak boson, while 
the photon can be emitted either from a fermion line or from a W boson in the i-channel. 
Representative Feynman diagrams for the production process are depicted in Fig. [2l The 
matrix elements for Hjjj production are combined with a generic decay of the Higgs boson 
into a fermion-antifermion pair. NLO-QCD corrections are supplied for the production process 
only. 

In order to speed up the calculation, we have used several approximations that have already 
successfully been applied to other WBF-type reactions. In our implementation, we disregard 
all diagrams of the type qq — >• VH'j with a subsequent V ^ qq decay, where V = Z, W. We 
also neglect interference terms between t- and n-channel diagrams. In the phase-space regions 
relevant for the Higgs search in the WBF channel, the disregarded contributions have been 
shown to be small [14j. 

2.2. NLO-QCD corrections 

The virtual corrections comprise the interference of one-loop diagrams with the Born amplitude. 
Due to color conservation, only selfenergy, vertex, and box corrections to either the upper or 
the lower quark line need to be considered. Within our approximations, contributions from 
diagrams with a gluon connecting the upper and the lower quark line do not interfere with the 
corresponding Born diagrams. 

The combination of all virtual corrections to quark line i gives rise to an amplitude of the 
form 

MV = Ms^C, (1^) V(l + .) (-| - ^ + c.,,,) + MV , (1) 

with Cp = 4/3 and Cvirt = vr^/S — 7 in dimensional reduction (cvirt = 7r^/3 — 8 in conventional 
dimensional regularization) . For each quark line i, the quantity Qi is related to the momentum 
transfer of the incoming to the outgoing fermion with momenta pin and pg, respectively, by 
Q| = -(Pin -Pfi)^. 

The finite parts A4y are evaluated numerically by means of a Passarino-Veltman tensor 
reduction |15j. For a small fraction of phase-space points, we encounter numerical instabilities 
which can be traced back to small values of the emerging Gram determinants. Such instabilities 
are monitored by testing Ward identities at every phase-space point. The Aiy contributions from 
phase-space points where these Ward identities are violated are removed. We have checked that 




(a) (6) 
Figure 3. Representative set of loop diagrams for H^jj production. 

tlie number of phase-space points that do not fulfill the Ward identities with a relative accuracy 
of at least 10~^ is less than 0.1 permille and that the numerical error due to disregarding these 
contributions is small. A detailed description of the calculation of the loop diagrams can be 
found in [16j. 

The real emission contributions are obtained from the LO diagrams by supplying all possible 
insertions of a gluon on the quark lines and considering crossing-related diagrams with a gluon 
in the initial state. We combine them with the virtual contributions by the dipole subtraction 
formalism of Catani and Seymour [17j . 

2.3. Checks 

Extensive checks have been made to verify our results. We have prepared three independent 
implementations of the entire calculation, which agree to better than one permille at the level of 
integrated cross sections. Furthermore, we have compared the leading order and real emission 
amplitudes with amplitudes generated by MadGraph [18' and found agreement to about 12 digits 
at a random set of phase-space points. The leading order cross section was shown to be in 
accord with the results of MadEvent [29J within the given accuracy of the program. A precise 
comparison to the results of Ref. [10] was difficult, as not all parameters of their calculation 
were explicitly listed in that publication. Nonetheless, we could reproduce their cross sections 
within a margin of 3 — 5%. 

In addition, we have checked that the real emission amplitudes are invariant under QCD 
gauge transformations. All contributions obey QED gauge invariance. We furthermore validated 
our implementation of the dipole subtraction formalism by testing that the dipole terms are 
approaching the real emission contributions both in the soft and collinear limits. 

3. Results 

For the results that we are presenting here, we are assuming a proton-proton collider with 
a center-of-mass energy of ^/s = 14 TeV. We are using the CTEQ6L1 parton distribution 
functions at LO and the CTEQ6M set with as{mz) = 0.118 at NLO [20]. As electroweak input 
parameters, 

mz = 91.188 GeV, 

mw = 80.398 GeV , (2) 
Gf = 1.166 X 10"^ GeV-^ 

are taken. Throughout, we have set the Higgs mass to mn = 120 GeV. Fermion masses are 
neglected and contributions with external top quarks are entirely disregarded. 



Table 1. Cross sections (in [fb]) within the WBF cuts of Eqs.Q-dS]), for different scale factors 
^ = £,R = S,F- The statistical errors of the results are below the permille level and therefore not 
given here explicitly. 
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We are using the kx algorithm as described in [21] with a cone parameter of Rc = 0.7 to 
define a jet. To avoid collinear singularities associated with the final-state photon, we utilize the 
isolation criterion introduced by Frixione [22j . An event is considered acceptable if the hadronic 
energy deposited in a cone around the photon is limited by 

El — cos AR , , « ^ ^ 

PTi < ^PT7 VAi? < 5o. 3 
1 — cos On 

i:AB.i^<AR " 

Here, the summation index i runs over all final-state partons found in a cone of size AR in the 
rapidity-azimuthal angle plane around the photon, pxi denotes the transverse momentum, and 
ARi^i the separation of parton i from the photon, while 5q stands for a fixed separation. 

While we account for the kinematic distributions of the Higgs decay into two generic massless 
fermions (denoted by h below, although not necessarily indicating a bottom quark), the branching 
ratio for this decay is not included in the results we are presenting here. 

Our inclusive set of cuts is given by 

PTi > 20 GeV , m]f > 100 GeV , 

\yj\<5, 1^7,61 < 2.5, (4) 

ARik>OA, 5o = 0.7, 

where pTi stands for the transverse momentum of any jet, Higgs decay particle, or photon, and 
ARik is the separation in the rapidity-azimuthal angle plane of any two objects in the recombined 
final state. 

Another, more stringent, set of cuts which will be entitled as "WBF cuts" in the following 
additionally comprises the typical cuts for Higgs boson searches in WBF. For this scenario, 
two tagging jets of large invariant mass have to be widely separated in rapidity in different 
hemispheres of the detector. The decay products of the Higgs boson are required to be located 
within the rapidity gap between the two tagging jets: 

Ai?ifc >0.7, 
m]f > 600 GeV , 

\yji - yj2\ > 4, (5) 
Vji ■ yj2 < , 
yr <y,,b<yr''. 

Numerical results for the integrated cross section at LO and NLO QCD within the latter set 
of cuts are shown in Tab. [TJ For the renormalization and factorization scales, 

/^iJ = fJ-o, IJ-F = (.F /"o, (6) 
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Figure 4. Scale dependence of the integrated cross section within the inclusive cuts of Eq. Q 
at LO and NLO for two different choices of ^q. Shown are curves for = at LO (black dots), 
^j^ = = ^ at NLO (red sohd), = ^, = I at NLO (blue dash-dot), and = ^, ^p = 1 
at NLO (green dashes). 



we employ two different central scales, defined by /Xq = and /Xq = mj^ + YlPrj^ where the 
sum extends over defined jets. The scale factors, (,r = (,f = ^, ai^e varied by factors of two. 
Obviously, the dependence of the cross section on the factorization scale goes down when NLO- 
QCD corrections are included. The renormalization scale enters for the first time at NLO, since 
no strong coupling occurs in the tree-level process. For every scale choice, the ratio of the LO 
to the NLO cross section, referred to as K-factor, deviates from one by only a few percent. 

Figure H] illustrates the dependence of the integrated cross section within our inclusive set 
of cuts, Eq. ([U, on the unphysical renormalization and factorization scales and fip- For 
fiQ = Qf and fi^ = mjj + "^Ptj the curves follow the same trend. In each case, the scale 
dependence decreases drastically when NLO corrections are taken into account. In the range of 
0.5 < = C-R ^ 2, the value of the NLO cross section varies by just 2% for = Qf, and by 
1% for /iQ = m'jj + J2PTj- Comparing the results of Fig. [J] with those of Tab. [TJ we find that 
the WBF cuts of Eq. ([5]) reduce the signal by about a factor of two. 

For the selection of WBF events in a collider environment, understanding the kinematic 
properties of the tagging jets is of vital importance. Fig. [5] (a) shows the distribution of 
the hardest tagging jet at LO and NLO within the WBF cuts of Eqs. (jlD-(I5]). While the 
NLO corrections are small in total, the shape of the transverse-momentum distribution changes 
significantly when going from LO to NLO, with the peak of the distribution being shifted to 
slightly smaller values. In order to quantify the impact of the NLO-QCD corrections on the 
distribution of an observable O together with the scale uncertainties of the LO and the NLO 
prediction, we introduce the relative variation 

da{0/dO 

^(0) = d.NLO(^ = l)/dO-l' (7) 

where da{£,)/dO denotes the LO or NLO expression, evaluated for arbitrary values of the scale 
parameter i = = ^R- The choice of /^o is identical for da/dO and da^^^/dO. For the 
transverse momentum distribution, this relative correction is displayed in Fig. [5] (b) for two 
different values of ^. The bands in between the curves for = 2 and ^ = 0.5 indicate the 




Figure 5. Transverse momentum distribution of the hardest tagging jet at LO (dashed black 
hue) and NLO (sohd red hue) [panel (a)] and relative corrections according to Eq. ([7|), when 
the factorization and renormalization scales are varied in the range Qi/2 < = < 2Qi 
[panel (b)]. The WBF cuts according to Eqs.dH and ([5]) are employed here. 
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Figure 6. Invariant mass distribution of the Higgs boson-plus-photon system at LO (dashed 
black line) and NLO (solid red line) [panel (a)] and relative corrections according to Eq. ([7]), 
when the factorization and renormalization scales are varied in the range Qi /2 < = /i_F < 2Qi 
[panel (b)]. 



scale uncertainties of the LO and the NLO predictions, respectively. At NLO, this uncertainty 
is smaller than 4% over the entire pT range considered. The scale dependence is much more 
pronounced at LO, though, and is generally larger for higher transverse momenta. 

In order to distinguish the Higgs signal in the H ^ bb decay mode from QCD-induced 66(7) 
backgrounds, imposing a cut on the invariant mass of the bbj system can help [lOJ. Figure O 
illustrates the perturbative uncertainties associated with the invariant mass distribution of the 



Higgs boson-plus-photon system, being reconstructed from the four-momenta of the photon and 
the decay products of the Higgs boson. The distribution vanishes for Mh-^ < mn and peaks 
for rriH = 120 GeV at around Mn-y ~ 165 GeV. For larger Higgs masses, da/dMn-y would be 
shifted to correspondingly higher values. 



4. Summary 

We have presented an NLO-QCD calculation for Higgs production in association with a 
photon via weak boson fusion at the LHC in the form of a fully flexible parton-level 
Monte Carlo program. Our calculation has been implemented in the publicly available 
VBFNLO package [T2l[T3]. 

Analyzing the pp — )■ H'jjj process for two different combinations of selection cuts, we 
have found that the impact of the NLO-QCD corrections on the integrated cross section is 
generally small. Its actual numerical value slightly depends on the choice of renormalization and 
factorization scales. The shape of some distributions changes, however, beyond leading order. 
For observables and regions of phase-space within typical WBF cuts, the NLO corrections are 
mostly smaller than 15%. 
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